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Johan Bijnens 



In this talk I will restrict myself to the light pseudo-scalar mesons (two, three or more flavours) 
and strong interaction only, including interaction with external currents and densities. This ex- 
cludes the work including baryons, heavy quarks, vector mesons, structure functions and related 
quantities and non-leptonic weak interactions as well as effects of photon loops. I start by men- 
tioning a few important historical papers which have roughly a jubilee this year. Section |l| followed 
by a short introduction to Chiral Perturbation Theory (ChPT) including discussions on chiral log- 
arithms and what exactly to expand in. The next two sections give an overview of two and three 
flavour ChPT respectively, where in the latter case I go into some detail into how the LECs are de- 
termined experimentally and the assumptions on the order LECs. Section ^ discusses the decay 
T] ^ 371 while at the end I make some comments about recent work in partially quenched ChPT 
and applications of the renormalization group to ChPT. Many part are also in my earlier talk 



1. Some History: 50, 40, 35, 30, 25, 20 and 15 years ago 

Since in this conference we celebrate the scientific achievements of Gerhard Ecker and Jiirg 
Gasser, it is appropriate to look back at the history of Chiral Perturbation Theory. I have picked 
out a few papers which feU at or close to jubileum years. About 50 years ago the subject was 
started with the Goldberger-Treiman relation [Q] and the advent of PC AC, the partially conserved 
axial-current [Q], and how this reproduced the Goldberger-Treiman relation. About 40 years ago 
a lot of work had been done within the framework of PCAC but 1968 and 1969 saw some very 
important papers: the Gell-Mann-Oakes-Renner relation [Q] and the proper way how to implement 
chiral symmetry in all generality in phenomenological Lagrangians Shortly afterwards loop 
calculations started with e.g. loop results for nn scattering and rj ^ 3n [^. 30 years ago 
the start with the modern way of including higher order Lagrangians and performing a consistent 
renormalization came with [^]. At the same time there was also the beautiful paper by Gasser and 
Zepeda about the types of non-analytical corrections that can appear [^. The seminal papers by 
Gasser and Leutwyler of 25 years ago then put the entire subject on a modern firm footing [ ]T0| , 



11]. The same period also had my own entry into the subject [12]. Lots of one-loop calculations 



were done and the understanding that the coefficients in the higher-order Lagrangians could be 



understood from the contributions of resonances was put on a firm footing [ ]10| , |13| ]. Let me close 
this historical part with two 15 year old papers, a very clear discussion of the basics of ChPT [14] 
and the first full two-loop calculation [[T^]. 



2. Chiral Perturbation Theory: ChPT, CHPT or ;t;PT 

ChPT is best described as " Exploring the consequences of the chiral symmetry of QCD and 
its spontaneous breaking using effective field theory techniques" and a particularly clear discussion 
about its derivation and underlying assumptions can be found in [14]. Some reviews are [ |l6| , |I7| ]. 



More reviews and references to introductory lectures can be found on the webpage [ |18[ ] 



For effective field theories, there are three principles that are needed and for ChPT they are 

• Degrees of freedom: Goldstone Bosons from the spontaneous chiral symmetry breakdown. 

• Power counting: This is what allows a systematic ordering of terms and is here essentially 
dimensional counting in momenta and masses. 
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• Expected breakdown scale: The scale of the not explicitly included physics, here reso- 
nances, so the scale is of order Mp , but this is channel dependent. 

Chiral symmetry is the (continuous) interchange of quarks. If we look at the QCD Lagrangian 
^QCD = £ [iqi^qi + m^qR - {qRqt + qiqR)] - ^Gl^G"^" (2. 1) 

q=u,d,s 

we see that we have an SU (3)\/ symmetry for equal quark masses but for = we can change left- 
and right-handed separately giving a SU (3)l x SU (3)r symmetry with q^ = {u d s) transforming 
as qt giqt and qR gRqR. 

The chiral symmetry is broken spontaneously by vacuum condensates (qq) = {qiqR + qRqt) 7^ 
0. This breaks the eight axial generators of the symmetry group but leaves the vector part unbroken: 
SU (3)l X SU {3)r—^ SU {3)v- This produces eight massless Goldstone Bosons and their interaction 
vanishes at zero momentum. The latter is very important, it is the reason why there exists a proper 
power-counting in ChPT. This is illustrated in Fig. 

p2 V Y {p')H^/p'?p'=p' 



{p^){Vp-)p' = p' 




i/p2 

Sd'p P' 




Figure 1: An illustration of the power-counting in ChPT. On the left we have the lowest order vertex with 
two powers of momenta or masses, the meson propagator with two inverse powers and the loop integration 
leading to four powers. On the right hand-side we see two one-loop contributions and how the counting on 
the left leads to the same power p"^ for both diagrams. This counting can be generalized to all orders 

We now start to look at the needed Lagrangians. The SU (3)l x SU {3)r/SU (3)y manifold is 
parameterized by a matrix 



U = exp (^/a/(2)<I>/Fo) with ^{x) 



(2.2) 



The traceless Hermitian matrix <I> is written in the usual pseudo-scalar fields. With the covariant 
derivative D^U = d^U — ir^U + iUl^ , which includes the left and right external currents: r{l)^ = 
+ {—)a^ and the matrix x = 2Bo{s + ip) that contains the external scalar and pseudo-scalar fields 
s and p, the lowest order Lagrangian is 

■^2 = {Fi/4) { {D^U'^D^'U) + ix^'U + XU'') } . (2.3) 

(A) is the trace over flavours Trp (A). Quark masses are included via s = diag(mi,,mt/,m^) H . 

At higher orders the number of terms in the Lagrangian increases rapidly. There exist two 
types, those representing contact terms, i.e. without pseudo-scalar bosons, and those with. The for- 
mer can never be measured but are the reflection in ChPT of the definition of currents and densities 
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2 flavour 3 flavour 3+3 PQChPT 

F,B 2 Fo,Bq 2 Fq,Bq 2 

/ 7+3 L':,H[ 10+2 L';,H[ 11+2 

/ c't 52+4 Cf 90+4 112+3 



Table 1: The number of parameters+contact terms for the various types of ChPT. 



in QCD. The latter are usually called low-energy-constants (LECs). The number of parameters at 
the various orders is shown in Tab. [I[ The order is from [19], order p'* from [ |l0| , |ll[], order p^ 



from [20 1 after an earlier partial result ^A\\. The partially quenched results are derived from the np 



flavour case [22]. The difficulty in obtaining a minimal set can be seen from the recent discovery 
of a new relation for two flavours [^]. Since the normal case is a continuous limit of the partially 
quenched case, the resulting LECs are linear combinations of partially quenched LECs using the 
Cay ley-Hamilton relations given in [^^. The general divergence structure at this order is known 
[^4|]. The parameters B ^ Bq and F 7^ Fq are the two versus three-flavour lowest order constants, 
these are different quantities. 

The main predictions of ChPT are twofold. 1) It relates processes with different numbers of 
pseudo-scalars. 2) It predicts nonanalytic dependences at higher orders, often referred to generi- 
cally as Chiral Log(arithm)s. As an example, the pion mass for = 2 is given at NLO by [ [lO| ] 



m 



I = 2Bm + 



1 , {2Bm) , 



3271^ ° n 

The implicit /i dependence in and the explicit dependence in the logarithm cancel. 



+ ••• (2.4) 



The LECs, like in ( |2.4j ), have to be determined experimentally or from lattice calculations. 
For np = 2 Ref. [[I^ introduced the /i independent li = {Ji2k^ /yi) /[(/x) — log (M^//x^) , which are 
proportional to the LECs Vj{iJL = mji). For np = 3 some of the corresponding yi are zero and no 
good equivalent definition of L,- exists. Here we always quote the L[(jLt). The scale jj. is arbitrary 
but becomes relevant when using estimates for higher order constants. 

A question which is often misunderstood is what quantities to expand in. The ChPT ex- 
pansion is in momenta and masses. However, one first has to decide whether to expand in low- 
est order quantities, like F,2Bm, or physical masses and decay constants, like mji,mK,m-q,Fji,FK- 
The latter is not unique either since the Gell-Mann-Okubo relation and kinematical relations like 
s + t + u = 2m\ + 2m\ for TT^-scattering can be (and are heavily) used to rewrite expressions. This 
sounds trivial but can change much how a series convergence looks as shown below for a simple 
example. I prefer to use physical masses and decay constants rather than the lowest order quanti- 
ties. The physical quantities are typically better known and the chiral logs are created by particles 
propagating with their physical momentum. Also, thresholds appear in the right places at each 
order in perturbation theory. The differences are higher order, but can be numerically important. 

Take as a simple example the relations m^ = mQ/ {\+ arriQ //o), '«;r = /o/(l + bniQ //o ) , as 
exact. We can expand to NNLO in several ways, some examples are 

mk -jirin ( mo -^mk \ 

m„ = mo-a-^+a^-^ + --- = fJ\ - +b^-^ + ■ ■ ■] (2.5) 

Jo /o V /o /o / 
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ntj, 



mo 



-a^ + a{b- 
Jn 



JTC 



fn = fo{l-b'^ + b{2b- 

In 



-a)-f + . 

J 71 



(2.6) 



The coefficients in the expansion and the actual numerical values clearly depend on the way we 
write the results. The plots in Fig. ^ show the convergence for a = I, b = 0.5 and /o = 1. Only 





0.2 0.3 
mo 



Figure 2: On the left the expansion of m^i in terms of m„//o (2.5) and on the right in terms of niji/ f,i (2.6) 
Shown are the full results (m^i) and the first three approximations. 



knowing the first three terms one would draw very different conclusions on the quality of the 
convergence from Fig. |2|for the different ways of expanding. 



3. Two-flavour ChPT at NNLO 



References to order and p'^ work can be found in [[l7|]. The first work at NNLO used 
dispersive methods to obtain the nonanalytic dependence on kinematical quantities, q^,s,t,u at 



NNLO. This was done for the vector (electromagnetic) and scalar form-factor of the pion in [ |25| ] 
(numerically) and 



(analytically) and for TTTi-scattering analytically in []27p. The work of [ |27p 
allowed to put many full NNLO calculations in two-flavour ChPT in a simple analytical form. 

Essentially all processes of interest are calculated to NNLO fully in ChPT starting with 77 — > 
71° 71° [|T5|, 77 ^ 71+ 71" l|2|, |33l, and Mj, 7r7r-scattering [Q, the pion scalar 



and vector form-factors [32] and pion radiative decay n Ivy [33]. The pion mass is known at 
order p^ in finite volume |34]. Recently 71° 77 has been done to this order as discussed in the 
talk by Kampf 



The LECs have been fitted in several processes. I4 from fitting to the pion scalar radius [ 33 , 36], 



Ij, from an estimate of the pion mass dependence on the quark masses [|10|, ^ and h , h from the 
agreement with TTTT-scattering [36], l(, from the pion charge radius [ |3^ ] and le — 15 from the axial 
form-factor in 71 — > Ivy. There is also a recent determination of I5 from hadronic tau decays [^]. 



The final best values are [|3|,|3|,|3|,|3] 



/"i = -0.4 ±0.6, 
Z6-/5 =3.0±0.3, 



/2= 4.3 ±0.1, 

k = 16.0 ±0.5 ±0.7, 



Z3 = 2.9±2.4, Z4 = 4.4±0.2, 
/5 = 12.24 ±0.21. 



(3.1) 



There is information on some combinations of p^ LECs. These are basically via the curvature 



in the vector and scalar form-factor of the pion [32] and two combinations from 7r7r-scattering [36] 
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LO 
NLO 
NNLO 



0.1 0.15 
[GeV^] 



0.15 h 

0.1 
0.05 




LO 
NLO 
NNLO ■ 



0.1 0.15 
[GeV^] 



Figure 3: The pion mass squared as a function of the quark mass viaM^ = 2Bm, left with inputs as in (3.1 ) 
and right with lo, — 0, both are for np —2 ChPT. 
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Figure 4: The pion decay constant as a function of the Figure 5: The tttt scattering length in three- 
quark mass via = IBin, for n/r = 2 ChPT. flavour ChPT as a function of the input values of 



L4,Lg used in the fits, from 



from the knowledge of and be in that reference. The order LECs c'j are estimated to have a 
small effect for m„,f„ and TTTT-scattering. 

Let me now show the dependence on the quark mass via = 2Bm for m| with surprisingly 



small NLO and NNLO corrections for the values of the input parameters in (3J_) and cj^(/x = 
0.77 GeV) = 0. The full result is extremely linear as can be seen in the left plot in Fig. ^. The 
linearity is a consequence of the fitting parameters as can be seen in the right figure in Fig. ^. 
Similarly, Fj^ as a function of expanded as in (2^) is shown in Fig. 0. The values of m|, F;i and 
are determined self consistently via an iterative method from the ChPT formulas quoted in [|^]. 



4. Three-flavour ChPT 



4.1 Calculations 

In this section I discuss several results at NNLO in mesonic three-flavour ChPT. The formu- 
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las here are much more involved than in two-flavour ChPT and while the expressions have been 
reduced to a series of well-defined two-loop integrals, the latter are evaluated numerically. Both 
are the consequence of the different masses present. The vector two-point functions [|^, ^] and 
the isospin breaking in the pco channel [^] were among the first calculated. The disconnected 



scalar two-point function relevant for bounds on L4 and Lg was worked out in [ |42| ] The remaining 



scalar two-point functions are known but unpublished []43p. Masses and decay constants as well as 
axial-vector two-point functions were the first calculations which required full two-loop integrals, 
done in the n and rj [39, 44] and the K channel [^]. Including isospin breaking contributions to 



masses and decay constants was done in [45|. After had also been evaluated to NNLO |46] a 
fit to the LECs was done as described below. The vacuum expectation values in the isospin limit 
were done in [46], with isospin breaking in [^] and at finite volume in [47]. 



Vector (electromagnetic) form-factors for pions and kaons were calculated in [ ]48| , [49| ] and in 
49|] a NNLO fit for Lg was performed. L[q can be had from hadronic tau decays or the axial 
form-factor in 7i,K ivy. The NNLO calculation is done, but no data fitting was performed [|50|]. 
A rather important calculation is the form-factor. This calculation was done by [ ]5T| , |52[ ] and 
a rather interesting relation between the value at zero, the slope and the curvature for the scalar 



form-factor obtained []51|]. Isospin-breaking has been included as well [[ 



Scalar form-factors including sigma terms and scalar radii [54] and %% [38] and Tr/T-scattering 



[55] have been performed as well and used to place limits on L\ and Lg. Finally, the relations 
between the l\^c\ and L\.,C\ have been extended to the accuracy needed to compare order results 



in two and three-flavour calculations [56] and there has been some progress towards fully analytical 
results for m| [57] and TiA'-scattering lengths [58]. The most recent results are T] ^ 371 [|59|], isospin 
breaking in Kc^ [|^. 



4.2 The fitting and results 



The inputs used for the fitting, as discussed more extensively in [H3, 46], are 



• KfA. F{0), G(0), A from E865 at BNL[]60|]. 

• m^o> nt^+, m^, electromagnetic corrections include the violation of Dashen's theorem. 

• andFK+/Fj,+ . 

• nis/tn = 24. Variations with rris/tn were studied in [45, ^]. 



• L'^,L'^ the main fit, 10, has them equal to zero, but see below and the arguments in []42|]. 



Some results of this fit are given in Tab. g. The errors are very correlated, see Fig. 6 in [ |46| ] 
for an example. Varying the values of L^4,L'^ as input can be done with a reasonable fitting chi- 
squared when varying 10^ L4 from —0.4 to 0.6 and Lg from —0.3 to 0.6 [p]]. The variation of many 
quantities with L^jLg (including the changes via the changed values of the other LJ") are shown in 
[54, 38, Fit B was one of the fits with a good fit to the pion scalar radius and fairly small 



corrections to the sigma terms [ 54 ] while fit D [^] is the one that gave agreement with nn and 
TT^-scattering threshold quantities. 

Note that niu/md = is never even close to the best fit and this remains true for the entire 
variation with Z.4,Lg. The value of Fq, the pion decay constant in the three-flavour chiral limit, can 
vary significantly, even though I beheve that fit B is an extreme case. 
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fit 10 


same 


fitB 


fitD 




0.43 ±0.12 


0.38 


0.44 


0.44 




0.73±0.12 


1.59 


0.60 


0.69 




-2.53 ±0.37 


-2.91 


-2.31 


-2.33 




= 


= 


= 0.5 


EE 0.2 




0.97±0.11 


1.46 


0.82 


0.88 




= 


= 


= 0.1 


= 




-0.31 ±0.14 


-0.49 


-0.26 


-0.28 




0.60±0.18 


1.00 


0.50 


0.54 




5.93 ±0.43 


7.0 


- 


- 


IBoih/ml 


0.736 


0.991 


1.129 


0.958 


2 4 6 


0.006,0.258 


0.009,= 


-0.138,0.009 


-0.091,0.133 


2 4 6 


0.007,0.306 


0.075,= 


-0.149,0.094 


-0.096,0.201 


2 4 6 


-0.052,0.318 


0.013,= 


-0.197,0.073 


-0.151,0.197 


mu/nid 


0.45±0.05 


0.52 


0.52 


0.50 


Fo [MeV] 


87.7 


81.1 


70.4 


80.4 


£k. „4 6 

f^ - P ,P 


0.169,0.051 


0.22,= 


0.153,0.067 


0.159,0.061 



Table 2: The fits of the LJ" and some results, see text for a detailed description. They are all quoted at 
H = 0.77 GeV. Table with values from [|| ||, g, |||l|]. 



In Fig. ^ I show how the threshold parameter depend on the variation with Z.4,Lg. always 
agrees well with the result of [36] while only agrees well within a limited region [|^. For 
comparison, the order values are Oq = 0.159 and = —0.0454. The planes in Fig. ^ indicate 



the results 



0.220 ±0.005, 



0.0444 ±0.0010 [g6p. The same study was performed for 
nK scattering lengths in []55[] with the results of the Roy-Steiner analysis []62[]. The resulting limits 



on the input values of ^4,^6 are shown in Fig. |^ The resulting region called fit D in Tab. ^ is 



TiTt constraints 



jiK constraints 





-0.4 -0.2 



0.2 0.4 0.6 



Figure 6: The bounds on L^^U^ from nn and Tr/T-scattering threshold parameters. Left nn where the bound 
from Oq shown in Fig. |^ is the most stringent. Right nK. White regions are allowed. The region of fit D, 
compatible with both, is indicated by the circle. From p^]. 

lO^L^ K, 0.2, lO^Lg PS 0.0. This general fitting obviously needs more work and systematic studies 
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0.2 0.4 0.6 0.8 1 

r"B'(r"B)phys 

Figure 8: For fit 10; Fk/Fj^ as a function of nis Figure 9: A schematic indication of the estimate of the 
with nis/rh fixed. order LECs by resonance exchange. 



and constraints from lattice QCD on ^4,^^ will be very useful. 

I now show the dependence of a few quantities on the input masses. These are updates of the 
plots shown in [|4^, more can be found in |[T]]. A selfconsistent set of m\, m\, m^, F^, Bonis and 
Bom with the fitted values of LJ' and Fq is determined for each input value of two masses. This is 
done by iterating the formulas till convergence is reached. I show m| for fit 10 and fit D keeping 
Ms/fn = 24 and varying in Fig. 0. The large corrections for fit 10 come from the kaon mass. 
The decay constants ratio Fk/Fji is shown as a function of with m,v/m = 24 as well. 



4.3 C]: estimates of order LECs 

Most numerical analysis at order use a (single) resonance approximation to the order p^ 
LECs. This is schematically shown in Fig. ^. The main underlying motivation is the large Nc 
limit and phenomenological success at order p^ [13]. There is a large volume of work on this, 
some references are [63]. The numerical work I will report has used a rather simple resonance 
Lagrangian 113, 31 , 45, 13] only. The estimates of the C\ is the weakest point in the numerical 
fitting at present, however, many results are not very sensitive to this. The main problem is that 
the C[ which contribute to the masses, are estimated to be zero except for r]' effects and how these 
might affect the determination of the others. The estimate is /i -independent while the C\ are not. 
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input 


cs-+4q 




q±3C3^ 


C[ + 4C'^ + 2C'2 




20.7 ±4.9 


-9.2 ±4.9 


9.9 ±2.5 


2.3 ±10.8 


'■ ^30'^ll 'Ql 


28.1 ±4.9 


-7.4±4.9 


21.0±2.5 


13.4±10.8 








23.5 ±2.3 


18.8±7.2 


Resonance model 


7.2 


-0.5 


10.0 


6.2 



Table 3: Different determinations of the same combinations of the C[ from nn and nK scattering [|65[]. 



The fits done here in 1 45, 46, 54] try to check this by varying the total resonance contribution 
by a factor of two, varying the scale jJ. from 550 to 1000 MeV and compare estimated C[ to exper- 
imentally determined ones. The latter works well, but the experimentally well determined ones are 
those with dependence on kinematic variables only, not ones relevant for quark-mass dependence. 

We are at present [ 64 ] working on a new fit and trying to find how can be done without these 
estimates. That there might be some strain can be seen from the different Q estimates from [65] 
shown in Table |3| using the results of 7171 and nK scattering of [38, 55]. 



5. rj 



KKK 



In the limit of conserved isospin, no electromagnetism and niu = m^i, the 77 is stable. Direct 
electromagnetic effects are small [^]. The decay thus proceeds mainly through the quark-mass 
difference mj, — mj. The lowest order was done in [67], order p'^ in [ |68| ] and recently the full order 
has been evaluated [59]. The momenta for the decay 77 —>■ n^n^n^ we label as Pn, p+, p- and 



po respectively and we introduce the kinematical Mandelstam variables s = + p-) ,t = {p+ + 
poY , u = {p- + po)-^ . These are linearly dependent, s + t + u = m\o + ± ± = 3^0 .The 



amplitudes for the charged, A{s^t,u), and neutral, A{s,t^u) are related 

A{si,S2,S3) = A{S\,S2,S^) +A{s2,Si,,Si) +A{si„S]_,S2) 



(5.1) 



The relation in ( |5.1| ) is only valid to first order in m„ — mj. The overall factor of m„ — can be 
put in different quantities, two common choices are 



A{s,t,u) = —M{s,t,u) or A{s,t,u) = -^-^{mj,-mj.) ^^^^ , 



(5.2) 



with R = {nis — m) jinid — niu) or = R{ms + nid) /(2m) pulled out. The lowest order result is 



M{s,t,u)Lo = [{^mml-s) /F^. 



(5.3) 



The tree level determination of 7? in terms of meson masses gives with ( |5.3| ) a decay rate of 66 eV 



which should be compared with the experimental results of 295±17 eV[|69|]. In principle, since the 
decay rate is proportional to 1//?^ or 1 / Q"^, this should allow for a precise determination of R and Q. 
However, the change required seems large. The order p^ calculation [^] increased the predicted 
decay rate to 150 eV albeit with a large error. About half of the enhancement in the amplitude 
came from TTTT rescattering and the other half from other effects like the chiral logarithms []68|]. The 
rescattering effects have been studied at higher orders using dispersive methods in and [J7T|]. 
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Both calculations found an enhancement in the decay rate to about 220 eV but differ in the way 
the Dalitz plot distributions look. This can be seen in Fig. |l^ where I show the real part of the 
amplitude as a function of s along the line s = u. The calculations use different formalisms but 
make similar approximations, they mainly differ in the determination of the subtraction constants. 



That discrepancy and the facts that in the dispersive estimate [ 73 1 was about half the full ChPT 

2.5 I , , , n , , n 2.5 I , , 1 ^ — H ' ' n 



1.5 



0.5 



-0.5 




P 

dispersive 



1.5 



1 



0.5 



-0.5 



s/m 




^ 



Figure 10: Left: Decay amplitude obtained by use of Khuri-Treiman equations [[70p along the line s = u. 
Right: Alternative dispersive analysis for the decay amplitude[[7l|]. Figs, from [Q, adapted from [7l| | 



calculation [46] and at order p"^ the dispersive effect was about half of the correction for 77 
makes it clear that also for this process a full order calculation was desirable. 
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Figure 11: Left: The amplitude M{s,t,u) along the line t — u. The vertical lines indicate the physical region. 
Shown are the real and imaginary parts with all parts summed up to the given order. Right: Similar plot but 
along the line s = u. Figs, from [p9[. 



The calculation [ p9| ] generalized the methods of | |45| ] to deal with n -rj mixing. Here I show 



only results. In Fig. [11| I show the numerical result for the amplitude along two lines in the Dalitz 
plot, t = u and s = u. The latter can be compared directly with the dispersive result of Fig. The 
correction found in [ ]59| ] at order p^ is 20-30% in amplitude, larger in magnitude than the dispersive 
estimates [ |70| , [7l| ] but with a shape similar to ||7l|]. 

The Dalitz plot in T] ^ 3 tt is parameterized in terms of x and y defined in terms of the kinetic 
energies of the pions 7] and 2rj = "Itj — 2mji+ — nij^o for the charged decay and z defined in terms 
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Exp. 


a 


b 


d 


KLOE 
Crystal Barrel 
Layter et al. 
Gormley et al 


-1.090±0.005+o:o?9 
-1.22 ±0.07 
-1.08±0.014 
-1.17±0.02 


0.124±0.006± 0.010 
0.22 ±0.11 
0.034 ±0.027 
0.21 ±0.03 


0.057±0.006lo:o?6 
0.06±0.04 (input) 
0.046 ±0.031 
0.06 ±0.04 



Table 4: Measurements of the Dalitz plot distributions in 7] ^ Tt'^n n^. Quoted in the order cited in 
The KLOE result / is / = 0. 14 ± 0.01 ± 0.02. 





^0 


a 


b 


d 


f 


LO 


120 


-1.039 


0.270 


0.000 


0.000 


NLO 


314 


-1.371 


0.452 


0.053 


0.027 


NLO {L] = 0) 


235 


-1.263 


0.407 


0.050 


0.015 


NNLO 


538 


-1.271 


0.394 


0.055 


0.025 


NNLO (jU = 0.6 GeV) 


543 


-1.300 


0.415 


0.055 


0.024 


NNLO ill = 0.9 GeV) 


548 


-1.241 


0.374 


0.054 


0.025 


NNLO iC] = 0) 


465 


-1.297 


0.404 


0.058 


0.032 


NNLO m = q = 0) 


251 


-1.241 


0.424 


0.050 


0.007 



Table 5: Theoretical estimate of the Dalitz plot distributions in T] ^ tt+tt tz . 



of the pion energies £",-. The amplitudes are expanded in x = \/3 (r+ — T^)/Q-^,y = ?>Tq/Q-^ — 1, 



Z = (2/3) Zi=\3 {^Ei -mr^Y / {mn - ^m^of, via 

\M{s,t,u)\^ = Al{l+ay + by^ + dx^ + fy^ + ---) , \M{s,t,u)\^ = Ag (1 +2a2 + • • •) -(5.4) 

Recent experimental results for these parameters are shown in Tabs. ^ and |^ There are discrep- 
ancies among the experiments but the latest precision measurements of a agree. The predictions 
from ChPT to order with the input parameters as described earlier are given in Tabs. ^ and ^. 
The different lines corresponds to variations on the input and the order of ChPT. The lines labeled 
NNLO are the central results. The agreement with experiment is not too good and clearly needs 
further study. Especially puzzling is that a is consistently positive while the dispersive calculations 
as well as [75] give a negative value. The inequality a < (^d-\-b — /4 derived in [59] shows 
that a has rather large cancellations inherent in its prediction and that the overestimate of Z? is a 
likely cause of the wrong sign for a. The fairly large correction gives in the end larger values of Q 
compared to those derived from the masses [B^ 



6. Even more flavours at NNLO (or PQChPT) 

NLO Partially Quenched ChPT has been studied by many people and found to be very useful, 



see [ l78[ ] and references therein. The masses and decay constants are known to NNLO for almost 
all possible mass combinations. Formulas were kept in terms of the quark-mass expansion to 
avoid the proliferation in physical masses appearing in this case. The three sea flavour masses and 
decay constants are in ^] and the two sea flavour results are in [pO|]. Numerical programs 
are available from the authors. The formulas are in the papers but can be downloaded from [jH 
PQChPT NNLO results for neutral masses are in 
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a 


Exp. 


a 




LO 


1090 


0.000 


Ciystal Ball (MAMI C) 


-0.032 ±0.003 




NLO 


2810 


0.013 


Crystal Ball (MAMI B) 


-0.032 ±0.002 ±0.002 




NLO m = 0) 


2100 


0.016 


WASA/COSY 


-0.027 ±0.008 ±0.005 




NNLO 


4790 


0.013 


KLOE 


-0.027 ±0.004;JJ;{J°^ 




NNLO (C[ = 0) 


4140 


0.011 


Crystal Ball (BNL) 


-0.031 ±0.004 




NNLO (L[,C[ = 0) 


2220 


0.016 


WASA/CELSIUS 


-0.026±0.010±0.010 




dispersive [[70|| 




(0.007—0.014) 


Crystal Barrel 


-0.052±0.017±0.010 




tree dispersive 




-0.0065 


GAMS2000 


-0.022 ±0.023 




absolute dispersive 




-0.007 


SND 


-0.010±0.021 ±0.010 




Borasoy Jt^ 




-0.031 






error 


160 


0.032 



Table 6: Measurements of the Dalitz plot distribution in 



J] — > Tt^n^n'^. Quoted in the order cited in |f7^. Table 7: Theoretical estimates of the Dalitz 

plot distribution in Tj ^ ;r'';r^;r*'. 1^ 



7. Renormalization Group 

In ChPT the renormalization group is not quite as useful as in renormalizable theories but 
Weinberg [||] already showed that some predictions can indeed be made. In particular leading 
logarithms at two loops can be had from only one-loop calculations. This was used for obtaining 
the leading double logarithms in nn scattering [ |82| ] and in general The extension to all orders 



was proven in []84[]. A leading log to five loops was obtained in [ ]85[ ] and recently recursion relations 
in the massless case were used to get results to very high orders [86]. 

The latter papers solved the practical problem of keeping track of all-order Lagrangians using 
two observations. First, in the massless limit, tadpoles vanish, and thus the number of external 
legs needed at any order does not increase. Second, the main vertex needed is then the four meson 
vertex. Here they found a useful general expressions using Legendre polynomials allowing them 
to do all needed loop integrals and obtain a fairly simple iterative algebraic recursion relation. 



8. Conclusions 



Modern ChPT is doing fine. Two flavour ChPT is in good shape: it is now precision science in 
many ways. For three flavour ChPT the corrections are larger and there seem to be some problems, 
but many parameters, especially in the scalar sector are rather uncertain and errors are very quantity 
dependent. Many partially quenched NNLO calculations have been done with an eye on lattice 
calculations and their extrapolations. A final comment is that new application areas continue to be 
found for ChPT and EFT in general. 
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